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Abstract 

Analytical periodic solutions for weakly Coupled Map Lattices are 
shown in an explicit form as well as in a recurrence relation. The re- 
sults establish a link between a matricial representation and recurrence 
relations of the solutions. 
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1 Introduction 



Regular higher order structures, chaotic or not, are very common in nature. 
They can be found in physics, chemistry, biology, social sciences and, in gen- 
eral, in any sytem with interacting elements. Regularity emerges as a con- 
sequence of the interaction among diverse elements of the system, therefore, 
to describe them properly, two things should be considered: the dynamics of 
every element in the system and the interactions among them. 

A common model used to simulate this situation is Coupled Map Lattices 
(CML), dynamical systems with discrete space and time, and continous state 
variables. Their origin is in the seminal work of Turing about morphogenesis 
[TJ. Later, Kaneko and coworkers studied them extensively (2} [3J SI [5J E]. 
CML have proved to be very good in describing collective behaviors observed 
in many fields: physics, chemistry, social sciences, biology [THE], that is why 
they are so extensively studied. On the other hand a continous non linear 
system can be studied by discretizing it [9], what yields a CML. A side effect 
of the discretization process is time-computer saving, that plays an important 
role when many equations are involved in collective phenomena. CML have 
also been used to study the control of spatiotemporal chaos [10J. 

Mathematically speaking, a one-dimensional CML [TT] is a 1-dimensional 
lattice where each site evolves in time according to 

X t {n + 1) = (1 - a)f{X t {n)) + s. ^ 
i = 1, m 

where Xi{n) denotes the state variable for the site i at time n. The parameter 
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a weights the coupling among sites. Here, periodic boundary conditions are 
asumed 

Xi(n) = X i+m {n) Vz 
For weak coupling, we write 

X^n + 1) = (1 - ea)/(Ai(n)) + f E?=i /(Xj (n)) 

Z = 1, • • ■ ,p £<1 

where a is the control parameter. 

In spite of the big amount of results in CML very few of them are an- 
alytical [12]. Therefore, patterns, periodic solutions, topological entropy or 
any other number that measures the complexity of a dynamical system can 
not be explicitly deduced. The goal of this paper is to find the analytical 
solution of (OQ), using perturbative methods and recurrence relations. The 
solution will be seeked in the form Xi{n) = x* + eAi where x* % — 1, • • • ,p 
are fixed points of f p and Ai is a perturbation term that we will explicitly 
calculate. 

2 Theoretical results 

In order to find the solution of (JTJ) we will use perturbative methods, because 
the system (pQ) is weakly coupled. So, we will seek a solution in the form 
Xi{n) = x* + eAi where f p {x*) = x* i = 1, • • • ,p. 

If we can find an expresion allowing us to know how the oscilator Xi(n) 
is turned into Xi{n + q) after g-iterations, then we will take advantage of the 



3 



periodicity of the solution to get A- We will firstly obtain this expression to 
reach our goal. 

Substituting x* + eAi into (PQ) gives 



ryf 

Xi { n + 1) = (1 - ea)f{x* + eAA + — £ f(x* + eA 3 ) 

V j= i 



(2) 



and expanding in powers of e (El) becomes 



X l (n+l) = f(x l )+e 



a 



AJ'(x*) - af(x*) + fix*) 



0(e 2 



or 



Xi(n + 1) = f(x*)+et h +0(e 2 ) 
t h = Afix*) - a f{x*) + a Y* =x f{x*) 



(3) 



By using ([3j) we obtain: 



X t in + 2) = (1 - ea)f(Xi(n + 1)) + f E?=i /(^-(n + 1)) 

= (1 - ea)f ifix*) + et n ) + f £? =1 f{f{x*)) + 0( £ 2 ) 

= Pix*) + eUj'ifix*)) - eapix*) + f Yg=i Pi*)) + 0(e 2 ) 

= fix*) + e (Uj'ifix*)) - apix*) + f Ef =1 fix*)) + 0(e 2 ) 

= fi x *)+et l2 +Oie 2 ) 



That is 



X,(n + 2) = / 2 (x*) + ^ i2 + 0(£ 2 



**> = /'(/«))* 



a/ 2 ta) + ?E?=i/ 2 (* 



(4) 



4 



The former result sugests that 

Xi(n + q) = P(x*) + et iq + 0(e 2 ) 
U q = f'(f q - l (x*))U q -i ~ aP(x*) + f EU f q ( x j) 

This will be proven in the following theorem. 
Theorem 1 The CML 

Xl ( n + 1) = (1 - ea)/(X f (n)) + f E?=i W) 
i = 1, ■ ■ ■ ,p e<l 

becomes 

X l (n + q) = r(x*)+et lq+ 0(e 2 ) (5) 
after (g — l)-iterations, where t iq satisfies the recurrence relation 



iW'to)-«/to) + ?E5=i/(a$) « =1 
*i, = S ( 6 ) 

/'(/ rl (^))Vi - + f E?=i 9 > 2 

Proof. We will prove the theorem by induction. 

i) By examination of ([3]) and PJ the validity of ((5j) and (J6]) is obtained 
for q = 1 and q = 2. 

ii) We assume the expression of if is valid for p = q and we will prove it 
is valid for p — q + 1. 
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Let us assume the induction hypothesis 



X i (n + q) = r(x*)+et iq + 0(e 2 ) 
Therefore, it results 

Xi(n + q + 1) = (1 - eamXiin + q)) + f E P j=1 + ?)) 

= (1 - ea)f(f«(x*) + st iq + 0(e 2 )) + f E P j= i f(f q (x*) + et iq + 0(e 2 )) 
= (1 - ea) (f(f q (x*)) + f'(f q (x*))et iq ) + f Ef =1 f(f q (x*)) + 0{e 2 ) 
= f q+ \xt) + f'(f q (x*))et tq - eaP+\xt) + f E P j=1 f q+ \x*) + 0(e 2 ) 
= P+\x*) + e (f'(f q (x*))et lq - aP+\x*) + f E?=i f q+ \x])) + 0(e 2 ) 
= p+i( x *)+et iq+1 +0(e*) 
as we wanted to prove. 

To get a complete solution of the problem we need an explicit (non- 
recursive) expression for t iq . This is obtained in the following theorem. 

Theorem 2. The solution of the recurrence relation 

a p 

K = HrtDlV. - + -E nx)) q > 2 

V j=1 

is given by 

k = a (nx*)) f - a (E q -=\ nxt) mil /vco)) + 

+f E q -=\ (eu nd, /'(/*(*»)) - *f q (xi) + f eLi f q (xt) 

Proof. We will prove the theorem by induction. 
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i) Let us see that the expression is valid for q = 2. 
By using ((6]) results 

U 2 = f'(f(x*))tn - «/ 2 (x*) + f EU f 2 ( x i) 

tn = fVfa)) (Atf'W-afte) + ^ = ,f{x*))-ap{x*)^Yf j=l f 2 {x*) 



P V 

U 2 = Arf\x*)fXf(x*))-af(x*)f(f(x*))+-f'(f(x*)) E f(x*)-af(x*)+- £ f(x*] 

P j= i P j= i 

(8) 

Subtituting q = 2 into (C]) gives 

t l2 = A(f 2 (x*))' - a(f(x*)f(f(xm + %(Y,Uf(x*)f'(f(x*))) ~ 
«/VD + fEf =1 / 2 (x?) 

This result coincides with the one obtained in (JBJ). 
ii) We assume the induction hypothesis and substituting into 

U q+1 = f(f q (x*))t iq - af«+\x*) + J EU f q+1 ( x j) 
results 
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%+l 



+ 



+ 



fk (x* 



/'(/*(*») [a (/*(**))' - a (ei=\ r «) urJ n /'(/' 

f'(f q (x*)) [f E*=i (E? = i /"(**) nLn /'(/*(<))) - «/"(^) + f Ef=i /'(a:,*) 

cm))' - « (/'(/*(*») el\ r«) mi / ; (/*(^))) 
f el\ (e?=i r(^) mi fv k te)j) - afWxfiwtt) 

= a (/^(x?))' - « (e*=i m) nLn /'(/*«))) 

+ f el\ (ej=i r (**) ntn /"(/*«))) - «/ , (/ 9 (^))/ 9 «) 



fc 



+ f ELi(E? =1 r(^)nL„ /'(/'«))) -a/* +1 (*.-) + f e? =1 / 9+1 (^) 



As we have previously indicated, now we can obtain by using theorems 
1 and 2 and the periodic boundary conditions. 

Theorem 3. Let ja^,^, ■■■■> x p} be a p-periodic orbit of a C 2 function / 
such that f p '{x*) ^ 1 i = 1, - ■ ■ ,p then the CML given by 



X^n +1) = (1 - ea)f(Xi(n)) + f E?=i /(*») 
i = 1, ■ • ■ ,p e < 1 



(9) 



has a p-periodic solution, given by 



+ j) = x* +j + £A i+j 
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where 



A = 



a 



1 " (f p (xt))' 



p-i 

E 

n=l 



i=l,---, P 

j = o,---,p-i 



V \ p-i 



1 * 



-x 



i+n 



+ -E^- II /'(***) + h*+-e 



p 



p 



;=1 



with periodic boundary conditions 



^i+p ~~ X i 



Vi 



Proof = x* + eAj is p-periodic orbit when Xj(n + p) = x* + eA 

By using theorem 2 it turns out that 

X l (n + p)= fP(xt)+e[Mf p M)' 

-a (El=\f n (xm P k =nf(f k (^))) 

+ f En=l (H=l /"(^) /'(/*«))) 

-a/*(aj) + f E f =1 /*(**,)] + C(£ 2 ) 

As 

the following equality is obtained: 
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x* + eAi 



x* + e 



-a (EiZ\x* +n nrJ n f'(f k (a 

^ p 2^n=l \2^j=l x j ilk=n J U \ x i 



-ax* + J ELi ^ 



Eliminating and separating terms of A, yields 



(/*(*?))' 



tt En=l X i+n Ilfc=n f'{f k { X i)) 
+ ^ Xm=l (Sj=l ^j) rifc=n f'( x i+k 

-ax* + J Ef=i x? 



(- 



+ p Z)j=l X j j Ilfc=n f'( X i+k 



-ax* + ™ Z)f=i x* 



and it is finally obtained that 



-4; 



a 



(f p ( x t))' 




' X i+n 



V 



3=1 



k=n 



-X, 



(10) 



The result is valid while 1 — (f p (x*))' is 0(1), so that (jBJ), obtained after 
expanding in powers of e, holds. The solution given by theorem 3 represents a 
propagating wave in the CML. Traveling waves in CML may explain ordered 
structures observed in nature [13J. 

The expresion (fTOl coincides with the one obtained in |17j (see theorem 
4 therein), thus providing an alternative way of facing CML. In [T7] a ma- 
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trix representation was given, meanwhile here the solution is obtained from 
a recurrence relation. Therefore a link between both representation is es- 
tablished. Two different approaches to the same problem provide a powerful 
framework to draw theoretical conclusions, in particular in areas like CML 
where most of the published results are numerical. 

The analytical approach allows to derive theoretical expressions for syn- 
chronizated states, travelling waves, period doubling cascades, intermittency, 
and so on [T7]. However, depending on the problem, recurrence relations or 
matricial representation may become a more or less useful tool in their anal- 
ysis. 

3 Discussion and conclusions 

We have found analytical periodic solutions for weakly Coupled Map Lattices, 
using recurrence relations. From recurrence relations an explicit solution 
is deduced, that coincides with the one obtained using matricial approach, 
therefore, the link between both techniques is established. The solution is 
very general because the result is valid for an arbitrary number of oscillators, 
and furthermore the individual dynamics of every oscillator is assumed to be 
ruled by an arbitrary C 2 function. 

Analytical reults are particularly important for two reasons: 

First.- Analytical expressions can obviously be managed with mathematical 
tools and new physical phenomena can be deduced from them. 

Second.- Spurious results, due to finite precision in computer simulations, 
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are avoided [MfT5l[T6]. 

Another couple of interesting facts must be pointed out. On the one hand, 
the recurrence relation allows for an exact and efficient numerical solution in 
a computer. On the other hand, as the solution is obtained for an arbitrary 
number of oscillators, limits tending to infinity can be calculated. This is 
an important point in the study of bifurcation cascades and of the onset the 
turbulence in fluids, where a finite number of elements is not adequate to 
described the phenomena. 
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